FISIST/19-2002/CFIF 
[h^h/0209192| 



Solar neutrino oscillations and bounds on neutrino 
magnetic moment and solar magnetic field 

E. Kh. AkhmedovS and Joao PulidJ 

Centro de Fisica das Interacgoes Fundamentals ( CFIF) 
Departamento de Fisica, Instituto Superior Tecnico 
Av. Rovisco Pais, P-lO^O-OOl Lisboa, Portugal 

Abstract 

If the observed deficit of solar neutrinos is due to neutrino oscillations, neutrino 
conversions caused by the interaction of their transition magnetic moments with the 
solar magnetic field (spin-flavour precession) can still be present at a subdominant 
level. In that case, the combined action of neutrino oscillations and spin-flavour pre- 
cession can lead to a small but observable flux of electron antineutrinos coming from 
the sun. Non-observation of these z/g's could set limits on neutrino transition moment 
H and the strength and coordinate dependence of the solar magnetic field B±. The 
sensitivity of the i^e fiux to the product fJ.B± is strongest in the case of the vacuum os- 
cillation (VO) solution of the solar neutrino problem; in the case of the LOW solution, 
it is weaker, and it is the weakest for the LMA solution. For different solutions, differ- 
ent characteristics of the solar magnetic field B±{r) are probed: for the VO solution, 
the v'e flux is determined by the integral of B±(r) over the solar convective zone, for 
LMA it is determined by the magnitude of B± in the neutrino production region, and 
for LOW it depends on the competition between this magnitude and the derivative of 
B I (r) at the surface of the sun. 
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1 Introduction 



The observed deficit of solar neutrinos [|I|] compared to tfie expectations based on tlie stan- 
dard solar model and the standard electroweak model is now firmly established to be 
due to non-standard neutrino properties. In particular, the SNO Collaboration has demon- 
strated 0] that a significant fraction of solar Vf, is converted into some other active neutrino 
species, which can be z/^, v^^ or Vr. The most plausible and widely accepted explanation 
of the observed solar neutrino deficit are neutrino oscillations; however, some alternative 
possibilities are not ruled out yet. One of them is neutrino spin-flavour precession ^ due 
to the interaction of neutrino transition (flavour off-diagonal) magnetic moments with the 
solar magnetic field. Unlike the ordinary neutrino spin precession 0, the spin-flavour pre- 
cession (SFP) can take place even if neutrinos are Majorana particles; in this case it converts 
left-handed v^. into right-handed z/^ or v^^ which would be in accord with the SNO findings. 
Neutrino SFP can be resonantly enhanced in matter 0, ||, very much similarly to the 
resonance amplification of neutrino oscillations, the MSW effect 0. 

SFP of solar neutrinos, both resonance and non-resonance, can very well account for 
the observed solar neutrino deficit. It yields an excellent fit of all currently available solar 
neutrino data (see, e.g., [10 - 18] for recent analyses), even somewhat better than that of the 
large mixing angle (LMA) oscillation solution, which is the best one among the oscillation 
solutions. However, to account for the solar neutrino data, SFP requires relatively large 
values of the neutrino transition magnetic moment, [i ~ lO^^^yU^ for peak values of the solar 
magnetic field -Bq ~ 100 kG. Although such values of /i are not experimentally excluded, 
they are hard to achieve in the simplest extensions of the standard electroweak model. 

In the present paper we shall be assuming that the neutrino transition magnetic moment 
and/or solar magnetic field strength are significantly below the values necessary for the SFP 
mechanism to account for the solar neutrino deficiency (though not completely negligible). 
Our assumption is that it is neutrino oscillations that solve the solar neutrino problem, 
while the SFP is present at a subdominant level. What can then the observable effects of 
SFP be? Its influence on the survival probability of solar z/g will be small and essentially 
indistinguishable from a small change of the neutrino oscillation parameters. However, the 
combined action of neutrino oscillations and SFP can lead to a qualitatively new effect 
which is absent when only oscillations or only SFP are operative - the production of a flux 
of electron antineutrinos |0, |19|, |2y, |2l|, Since all the currently favoured oscillation 

solutions of the solar neutrino problem - LMA and LOW MSW solutions and vacuum 
oscillations (VO) - require the solar neutrino oscillations to be driven by a large mixing 
angle [Q, an observable flux of solar v^, can in principle be produced. 

In the present paper we address the question of what can be learned about the neutrino 
transition magnetic moments /i and the solar magnetic field by studying the solar z/g's. In 
particular, we discuss the bounds on \i and the strength and coordinate dependence of the 
solar magnetic field that can be derived from the current upper limits on the solar v^, flux 
as well as from future experiments in case the flux of Ve from the sun is not observed. 
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Experimentally, z/g's have a very clear signature and can be easily distinguished from 
the other neutrino species. The main problem with detecting z/g's from the sun is the 
background of electron antineutrinos from nuclear reactors. This background is a steeply 
decreasing function of neutrino energy; it becomes negligible for E > {5 -8) MeV. Therefore 
only the solar neutrinos can contribute to an observable flux of solar z/g's and the energy 
interval to be studied is ^ (5 - 15) MeV. 

Neutrino magnetic moments can also manifest themselves through the additional con- 



tribution to the z/e scattering cross section in the solar neutrino detectors (see, e.g., |25| for 



recent discussions). However, these contributions can only be noticeable if > 10 ^^fiB- 



While such large values of fx are consistent with the current laboratory upper bounds [|6 



they exceed the astrophysical bounds < (1 - 3) xlO^^^yU^ ^7] by more than an order of 
magnitude. In our study we shall be assuming the astrophysical upper bounds to be satisfied 
and so shall neglect the effects of neutrino magnetic moment on neutrino detection. 



2 Probability of z/g production 

We shall assume neutrinos to be Majorana particles and consider their evolution under the 
combined action of oscillations and SFP. For simplicity, we shall discuss the case of just two 
neutrino flavours, z/g and z/^, and their antineutrinos. 

There are essentially two ways in which u^s can be produced: (1) the originally pro- 
duced solar z/g first oscillate into z/^, which are then converted into z/g by SFP; (2) solar 
Vg first undergo SFP and get converted into z/^, which then oscillate into z/g. This can be 
schematically shown as 

osc. SFP _ /I x 

J^eL ^ l^fML > l^eR , (1) 

SFP _ osc. _ /„\ 
I^eL > Z^Mi? ^ t^eR ■ (2) 

The oscillations and SFP in these two chains of conversions can either take place in the 
same spatial region, or be spatially separated. In the former case, the amplitudes of the 



processes (1) and (2) interfere. It was shown in [19| that the interference is destructive, 
leading to a significant suppression of the solar z/g flux even if the probability of SFP is 
large. The reason for this is CPT invariance from which it follows that the matrix of the 
Majorana-type transition magnetic moments is antisymmetric; this, in turn, implies that the 
amplitudes (1) and (2) are of opposite sign ||2^. The cancellation between the amplitudes 



(1) and (2) is exact when the corresponding intermediate states (z/^ and z/^) are degenerate, 
e.g., in vacuum; the degeneracy is lifted by matter and/or twisting magnetic field (i.e. by a 
magnetic field whose direction in the plane transverse to the neutrino momentum changes 
along the neutrino path) ||2^ , so inside the sun the z/g production is not completely blocked. 
Still, it is strongly suppressed 0, and so we shall concentrate on the Ue production in the 



^Except, possibly, for the LOW solution of the solar neutrino problem, see discussion in Sec. 3.2. 
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sequence of two spatially separated processes corresponding to eq. (2) pT|: first, SFP inside 
the sun converts solar v^l into v^r which then oscillate into Vf-R in vacuum on their way 
to the earth. The alternative possibility, corresponding to eq. (1), can be disregarded as 
the magnetic field in the region between the sun and the earth is negligibly small. The 
probability that a Vei born inside the sun will reach the earth as VeR is then 

P(z/eL l^eR) = P{l^eL ^ l^^lR ] Rq) " P{l^^lR ^ l^eR ] Res) 

^^esj , (3) 

where Rq is the solar radius, R^s is the distance between the sun and the earth, Am^ and 
6 are the neutrino mass squared difference and mixing angle, and E is neutrino energy. We 
shall now concentrate on the calculation of the SFP probability P{veL ^ ^fiR ] Rq)- 

The evolution of the neutrino system under the consideration is described by the fol- 
lowing system of equations 0, |22[ : 

^ = (Ye - C2S) UeL + S26 V^L + /i^^e*'^ V^R , (4) 

'eR = + 026) VeR " f^B^C'"'' Z/^L + 826 V^R , (5) 

/ 

I 

\iR 



iv'. 



^ ^'fiL = S26 VeL - ftBj_e"'' VeR + (V/, + 02(5) I/^L , (6) 
^ ^'uR = fJ'Bj_e~"'' VeL + S26 UeR - (Vf, - C2S) V^R . (7) 



Here K = \f2GF[Ne{r) - Nn{r)/2] and = y2G^[-iV„(r)/2] are the matter-induced 
potentials of electron and muon neutrinos, respectively, N^, and being the electron and 
neutron number densities; S2 = sin 26', C2 = cos 26', 6 = Am'^/AE. In eqs. (4)-(7) the 
angle defines the direction of the magnetic field B±{r) in the plane orthogonal to the 
neutrino momentum, B±{r) = |B_L(r)|. In the case of twisting magnetic fields (0 = 0(r)), 



the probability of UeR production inside the sun can be enhanced [^, since little is 
known about the possible twist of the solar magnetic field, we set = hereafter. This 
approximation is expected to reduce the calculated flux of solar Uers and therefore make 
our upper bounds on fiB± more conservative. 

Since we are interested in the situation when the SFP probability is small, it can be 
calculated in perturbation theory. To find the u^r amplitude in leading order in fiB±, 
one needs to use in eq. (7) the "unperturbed" by SFP value of the amplitude u^l- This 
means that the fiB± terms in eqs. (4) and (6) should be omitted. These two equations 
then decouple from the rest of the system and reduce to the standard two-flavour evolution 
equations describing neutrino oscillations in matter. Next, we take into account that, due to 
the partial cancellation between the amplitudes of transitions (1) and (2) discussed above, 
one can neglect the UeR amplitude inside the sun; we therefore omit eq. (5) as well as the 
UeR term in eq. (7). From eq. (7) one then obtains, up to an irrelevant phase factor, the 
following expression for the amplitude of the u^l —>■ v^r transition: 

AiueL - y,R Rq) = r\B^{r) v,^{r) 1^^-^^')'^'^' dr . (8) 
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Here is the coordinate of the point at which the VeL was produced in the sun, and the 
amphtude Vf.i{r) is determined solely by neutrino oscillations, i.e. has to be found from 
eqs. (4) and (6) with /ii?_L = 0. 

It is easy to check that for the LMA and LOW solutions of the solar neutrino problem the 
solar UeL ^ ^fiL oscillations are adiabatic; one therefore can use the adiabatic approximation, 
which yields 



I'eiir) = cos 9{ri) cos 9{r)e ^■I'ri^^'^^ + sin 6'(rj) sin ^(r)e ^•^'■i^^'^^ . (g) 



Here 



E,,2ir) = Tco, uj=J( ^^-^ - c,6y + {.3,6)^ , (10) 



and 6{r) is the mixing angle in matter at a point r, defined through 



cos2e(r) = ^ . (11) 

u 

Substituting (||) into (||) one finds 

A{u^L i^„R Re) = r'^f^B^ir) [cos^(ri) cose{r)e-'<^'^''^ + sin0(ri) sme{r)e-'3'^''^] dr , 

J Ti 

(12) 

where 

g,,2{r) = £{^^^^'C,5Tuj)dr'. (13) 

This expression is relevant for the LMA and LOW solutions of the solar neutrino problem. 

In the case of the VO solution, neutrino oscillations inside the sun are essentially blocked, 
so that instead of eq. (^ one finds from eq. (4) 

/ N -i {Ve-C25) dr' , ^ 

^eL{r) = e ■'^^ . (14) 



Eq. (pp then yields 



A{veL ^ y^R Rq) = ^^B^{r)e ■'n dr. (15) 



3 Calculations of the expected VeR flux 

Since only the solar neutrinos with the energies E > b MeV can contribute to an 
observable flux of solar antineutrinos, the expected flux of z/gi^'s in this energy region is 
<l>p^(i?) = ^SB{E)P{ueL —>■ i^eR,E). We shall now concentrate on the calculation of the 
transition probability PifeL — > t'efl, E) for the LMA, LOW and VO solutions of the solar 
neutrino problem. 
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In eqs. (|T2D and (|T^), the pre-exponential factors in the integrands are in general smooth 
functions of r, but the complex exponential factors are very rapidly oscillating functions. 
The oscillations are especially fast for the LMA solution because it corresponds to the largest 
values of Am^. The integrals of this type are notoriously difficult to calculate numerically 
- one needs very fine integration steps in order for the integrals to converge. For example, 
for the LMA step ~ 10 ^Rq or smaller is necessary. 

Fortunately, there are well developed and accurate approximate analytic methods of 
calculating such integrals (see, e.g., [^). For our purposes, analytic expressions also have 
an advantage of allowing one to directly relate the expected fiux of solar u^r to simple 
characteristics of the solar magnetic field. The integrals in eqs. ([T^) and ([l5|) are of general 
type 

1= f /(x)e-^»(^)rfx, (16) 



where f{x) is a smooth function of x and |5''(x)| is large except possibly in the vicinity of 
a finite number of points in the interval (a, b). The integrals of rapidly oscillating functions 
are in general strongly suppressed because the contributions of neighbouring points tend to 
cancel each other. The exceptions are the endpoints of the integration interval, for which 
there are no neighbouring points on one of the sides, and the points where g'{x) = 0, which 
correspond to the extrema of the phase g{x). In the vicinity of these (stationary phase) 
points the phase changes slowly and the corresponding contributions to the integral are not 
suppressed. These contributions can be found in the stationary phase approximation, which 
is the complex version of the steepest descent approximation. 

The contributions of the stationary phase points to the integrals of the type ([I6D are 



0{1/ ^J\g"\), whereas the endpoint contributions are in general 0{l/\g'\). It is easy to check 
that for the integrals in eqs. (|l^) and (^) and for the values of neutrino parameters relevant 
for the LMA, LOW and VO solution of the solar neutrino problem the condition 

\g"{x)\/g'{xf<^l (17) 

is always satisfied. Therefore, if there are stationary phase points in the integration interval, 
they are in general expected to give the dominant contributions to the integrals. However, as 
we shall show now, this does not happen in the cases of interest to us because the stationary 
points either do not exist or their contributions are strongly suppressed by nearly vanishing 
pre-exponential factors in (|16D . 

Let us check the stationary phase conditions for various solutions of the solar neutrino 
problem. If one disregards neutrino oscillations and considers only SFP, the amplitude of 
the UeL ^iJ,R transitions is given by eq. ([T5|). The stationary phase condition g' = then 
reduces to 

Ve + V^- 2C2S = . (18) 

This is nothing but the resonance condition for neutrino SFP in the small 6 limit 0, § . As 
was discussed above, neutrino oscillations inside the sun can only be neglected in the case 
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of the VO solution. Since for this solution the parameter 5 is very small (~ 10"^^ eV), the 
resonance condition (TH) is satisfied essentially at the surface of the sun, where the magnetic 
field strength is known to be very small (~ 10 - 100 G); therefore the stationary phase point 
plays no role in this case. 

For the LMA and LOW solutions, the stationary phase conditions g'12 = ^ reduce to 



V; + V„ - 2co6 



2V„. 



(19) 



which is the SFP resonance condition in the presence of mixing. It is easy to see that 
this condition cannot be satisfied for large values of mixing angles 6. Indeed, eq. (|19|) has 
solutions only when 

1 - V. 

(20) 



1 — Y 



where Ye is the number of electrons per nucleon in matter. Its value varies between 0.667 
and 0.868 inside the sun, so that condition (U^ requires sin^ 29 < 0.5, too small for both 
LMA and LOW solutions. Thus, there are no stationary phase points contributions to the 
amplitude of the UeL ^^iR transition in these cases as well, and one has to consider the 
endpoint contributions. 

To obtain the contributions of the endpoints of the integration interval to an integral of 
the type (|I6|) we integrate it by parts. Integrating two times one finds 



/(x)e-^3(^) dx 



J{x) , fix) f{x)g"{x) 



+ 



-ig{x) 



g\x) 



g'{x) 



+ 0{l/g\x) 



(21) 



From ([T7|) it follows that the third term in the brackets can always be neglected. We shall 
now consider the amplitudes of the v^l i^eR transition for all three solutions of the solar 
neutrino problem of interest. 



3.1 LMA solution 

The current best fit values of the neutrino parameters for this solution are Am^ ~ 6 x 10~^ 
eV^, sin26' ~ 0.9 We will be interested in neutrino energies £^ ~ (5 - 15) MeV, 

so that 5 ^ (1 - 3)xl0~^^ eV. It is easy to check that the contribution of the first term 
in the square brackets in eq. (|T2|) is always at least one order of magnitude smaller than 
that of the second term, and so we shall neglect it. We now apply formula (pl|). One has 
/' oc {d / dr) [sin 6 B±\ ~ sin 6B'j_ since in the adiabatic regime 6' is small. The comparison 
of the first and the second terms in the brackets in eq. (|21|) then shows that the first term 
dominates when 

B±{r^) » |5l(r.)|M(r.) ^ 10-'^i?o \B'^{n)\ . (22) 

We consider only the contribution of the neutrino production point rj since the magnetic 
field at the final point of evolution r = Rq is negligible. Let us introduce the scale height 
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for the solar magnetic field strength Lb = \B2^{dBj_{r)/dr)\~^, which is a characteristic 
distance over which the magnetic field varies significantly. Condition (P^ then can be 
written as 

LB(ri) > lO-^i?0. (23) 

Magnetic fields of scale heights as small as < lO~^i?0 can only exist over very short 
distances and so cannot lead to any sizeable SFP. We shall therefore consider large-scale 
solar magnetic fields which satisfy (p3|). 



0.005 




0.001 - 



I ' ' ' ' ' ' ' ' ' 1 

5 7 9 11 13 15 

E (MeV) 

Figure 1: Probabilities P{i'eL ^eR) corresponding to the v^l i^fiR amplitudes of eqs. (|l^) 
(wiggly curve) and (smooth curve) for the LMA solution. Magnetic field linearly decreasing 
from Bq = 5 X 10^ G at r = O.O5i?0 to zero at r = Rq and fi = W~^^hb were chosen. 

The maximum of production of solar neutrinos which are of interest to us corresponds 
to To ^ O.OSi?©; from eqs. (^) and (|T2|) we then obtain 

(24) 

ri=OmRQ 

where we once again omitted an irrelevant phase factor. In the energy region of interest, the 
corresponding probability varies by less than 20% (see fig. 1), and so in first approximation 



sin 6{ri)fiB_i 
9'2{ri) 
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we can replace it by its mean value. From eqs. (j^), (|2^ ) and (|T3D one then finds 

2 



Pii^eL i^eR) ^ 1-8 X 10"^° sin^ 20 



10 kG 



(25) 



where we have taken into account that the z/^r — u^r oscillations in the space between the 
sun and the earth are in the averaging regime. Eq. ( P^D is our final result for the LMA case. 

In fig. 1. the probabilities Pivei — t'ei?) corresponding to the Vei — > ^-'^J.R amplitudes 
of eqs. ([T^ ) and (|2^) are shown (the wiggly and smooth curves, respectively). The fast 
oscillations of the former are due to the interference between the two terms in eq. ([T^ , 
whereas the latter is smooth because in obtaining it we neglected the (subleading) first term 
in (p!2|). Notice that the oscillations described by the wiggly curve are in fact unobservable 
since they average out when one integrates over the neutrino production region or takes into 
account finite energy resolution of neutrino detectors. 



3.2 LOW solution 

The best fit values of the neutrino parameters for this solution are Am? ~ 10~^ eV^ 



sin 20 ~ 0.98 ||2J]. For the interval of neutrino energies of interest, the parameter 5 is in 
the range 5 ^ (1.7 - 5)xl0^^^ eV. Once again, the contribution of the first term in the 
square brackets in eq. (|1^) is much smaller than that of the second term, and so we neglect 
it. The analysis is similar to that in the LMA case, but there is one important difference. 
In the LMA case, the value of g2 changes only by about a factor of two from the neutrino 
production point to the surface of the sun (^'2-^0 varies between ~ 10^ and ~ 5 x 10^). 
In contrast to this, in the LOW case it changes by a large factor - from 5^2-^0 ~ 10^ at 
r ~ O.O5i?0 to (72-^0 ~ 10 at r = Rq. This difference in the behaviour of g'2 is a consequence 
of the fact that the values of 5 are much larger and so contribute significantly to (72 ^it all 
values of r in the LMA case whereas in the LOW case they dominate near the surface of the 
sun but are negligible in the solar core; the values of g2 at small r are mainly determined by 
the matter-induced potentials. Indeed, in the LOW case one finds from eqs. (|^) and (|10D 

g2\r=RQ ~ (1 - C2) 5, 5'2|r-=O.O5i?0 ^ (Ve + V^,) |r=0.05Ro • (26) 

The smallness of g'2 bX r = Rq implies that the integral in eq. (|12|) is now dominated by the 
surface of the sun. Since the magnetic field nearly vanishes there, one can expect the main 
contribution to the integral to come from the second term in brackets in eq. (0), i.e. from 
the term proportional to the derivative of the solar magnetic field rather than to the field 
itself. Indeed, the condition for the domination of the term oc B'j_ is 

B^{RQ)'^[\B'^{r)\/g'2ir)]\r=R^, or Lb{Rq) <^ 0.1 R^ . (27) 

Since the magnetic field at the surface of the sun is very weak, this condition is likely to be 
satisfied. For example, it is satisfied for any magnetic field profile which near the surface 
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of the sun has the form B± = Bo[l — (r/i?©)"] + Bi, provided that n > 10{Bi/ Bq); for the 
magnetic field at the solar surface Bi = 100 G and the peak value of the field Bq = 10 kG 
this requires ra > 0.1. In what follows we shall be assuming that condition ( ^Tf ) is satisfied. 
The amplitude of the UgL — > i^fiR transition is then approximately given by 

A{ueL - Re) ^ sin^(r,) sin^ ^'f^^'^^llf^ ■ (28) 

{I — C2) 

We shall parameterize the derivative of the magnetic field strength at r = Rq as 

with Bq the peak value of the field of the convective zone (O.T-R© < r < i?©), and k a 
parameter; n = 1 would correspond to, e.g., a linearly decreasing towards the surface of the 
sun field with the maximum at the center of the convective zone. The probability of the 
Vf,L — ^ I'eR transition is then 



, cos^e 2 / 3 X 10-1^ eV^ 

P{l^eL ^ l^eR) ^ 1.85 X 10"^ — ^ 

sm o \ 



H Bo 



10-i2^B 10 kG 



(30) 



where we have taken into account that the Ufj_ii —>■ VeR oscillations in the space between the 
sun and the earth are in the averaging regime, and that the mixing angle in matter at the 
neutrino production point 9{ri) ~ 7r/2. As an example, for the best fit values of neutrino 
parameters for the LOW solution, E = 10 MeV, = 10"^^//^, Bq = 100 kG and k = 1, 
eq. (ID yields P{ueL ^ i^cr) ^ 1-4 x lO^^. 



It should be noted that, while the probability P^Vei ^cr) based on the perturbation- 
theoretic VeL ^^iR amplitude of eq. (|12]) (fig. 2, solid curve) is very accurate, the simplified 
expression (|^) (dashed curve in fig. 2) is only correct within a factor of three or four. The 
reason for that is the following. In the case of the LOW solution the quantity (72-^© changes 
only by about a factor of two throughout the solar convective zone, being ~ 10 - 20, i.e. not 
too large. Therefore the contribution of the integration endpoint at the solar surface to the 
integral in (|12D is not much bigger than that of the convective zone, and the approximation 
( pop is only an order of magnitude estimate. 

In deriving eqs. ( pS] ) and ( pO] ) we have assumed, in addition to (|27|), that the contribution 
of the term proportional to the magnetic field at the neutrino production point can be 
neglected compared to the contribution of the term proportional to the derivative of the 
field at the solar surface. The condition for this is 

g±(r.) < ^"f'^^i^; \B'Ar=RQ. or B^{n)<l()'BoK. (31) 
(1 - C2Yd^ 

If it is not satisfied, the probability P^feL — ^ t'eii) in the LOW case will be approximately 
given by eq. (pSj), as it does in the LMA case. 
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We pointed out earlier that the production of UeR inside the sun is suppressed because of 
the partial cancellation between the amplitudes of the channels (1) and (2). This cancellation 
is, however, partly compensated in the case of the LOW solution by an enhancement due to 
the fact that is small and so is Ve in the solar convective zone, which is the region where 
the SFP mainly occurs. Therefore in the case of the LOW solution there is an additional 
channel of the v^r production - direct production inside the sun through neutrino oscillations 
and SFP. At the same time, the flux of v^rs generated via the mechanism that we discussed 
so far is somewhat reduced if one takes the direct production into account; as a result, the 
total VeR flux changes only slightly. We have checked by a numerical integration of the 
system of the differential equations (l)-(4) that our approximation, in which we neglect the 
direct production but consider unsuppressed production through the chain of the processes 
(2) yields a very accurate prediction of the total solar v^r flux. 



3.3 VO solution 

The best fit values of the neutrino parameters for this solution are Aim? ~ 4.5 x 10^^'' eV^, 
sin 2^ ~ 0.93 ||2^. The VeL t^fiR transition amplitude is given by eq. (pISl). For the interval 



of neutrino energies of interest, the parameter 6 is in the range S ~ (0.77 — 2.3) x 10"^'^ eV, 
i.e. is negligible compared to Vg + essentially everywhere in the sun. One can therefore 
put 5 = in eq. (pISD , i.e. the u^l i^^r transition probability is practically energy 
independent. 

In the case of the VO solution, the parameter g' that plays a key role for integrals of 
rapidly oscillating functions of the type (|16D is simply g' = Ve + V^. The quantity g'{r)RQ 
changes from a value ~ 10^ in the center of the sun to nearly zero at its surface, being < 1 
inside the convective zone. Thus, it is not legitimate to use the approximate expressions of 
the type (^) which require g'RQ to be very large everywhere in the integration interval. 
We shall therefore estimate the transition amplitude differently in this case. 

Let us first notice that the phase g{r) that corresponds to eq. ([15|) , 

9{r)= r{K + V,)dr', (32) 







first grows rapidly with r and then saturates and slowly approaches its asymptotic value 
because of the steep decrease of Ve + V^. The contribution to the integral ( P^ from the 



region where g{r) rapidly grows is strongly suppressed, and the main contribution comes 
from the region where the phase changes little. We can therefore adopt an approximation 
of retaining only the contribution of the region r > tq where the change of the phase Ag < 1 
and neglecting the phase change in this interval. Since it is the change of the phase and 
not its absolute value that matters for the transition probability, this amounts to merely 
replacing in this interval the integrand of eq. ([15|) by fiB±{r). The lower boundary of the 
new integration interval tq is defined from the condition 

^(i?0)-(7(ro)~l. (33) 
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As a result, we obtain 



Rq 

/i / B±(r) dr 



sm2 2^sm2 (^Re 



(34) 



Solving eq. (^) for ro we find tq — O.SlTi?©. We compared the VeL i^eR. probabilities 
obtained using this approximation with those based on the full Uei — ^ t'/^R amplitude ( [T5| ) 
for a number of magnetic field profiles and found that the accuracy of the approximation 
(|3|) is typically about 20%. 

Let us introduce the average magnetic field strength in the interval < r < Rq through 

Bi=- / Bi(r)dr, (35) 



. B^{r)dr, 



where Rq — tq = O.183i?0. Then eq. 



can be rewritten as 

1 2 



/i 



B 



10-12;XB 10 kG 



sin2 2^sin2 (^R, 
\ AE 



(36) 



4 Discussion 



We calculated the probability of production of solar z/g/j's assuming that the solar neutrino 
deficit is due to neutrino oscillations while the spin-fiavour precession caused by the inter- 
action of neutrino transition magnetic moments with the solar magnetic field is present as a 
subdominant process. We considered the SFP in perturbation theory and obtained analytic 
expressions for the transition probability P^Uei —>■ ^br) valid for the LMA, LOW and VO 
solutions of the solar neutrino problem. We compared these analytical expressions with the 
results of numerical integration of the system of differential equations (l)-(4) and found very 
good agreement in all the cases. 

For each of the solutions of the solar neutrino problem we then obtained simplified 
approximate expressions for P{veL ^eR)-, which allowed us to relate this probability with 
simple characteristics of the solar magnetic field Bj^{r). For different solutions, different 
characteristics of the solar magnetic field -B±(r) are probed: for the VO solution, the Ve fiux 
is determined by the integral of B±{r) over the upper 2/3 of the solar convective zone, for 
LMA it is determined by the magnitude of B± in the neutrino production region, and for 
the LOW solution it depends on the competition between this magnitude and the derivative 
of Bx^{r) at the surface of the sun. 

The accuracy of the simplified expressions for P{veL J^ei?) is also different for different 
solutions: the error is less than 3% for the LMA solution and about 20% for the VO solution, 
while for the LOW solution the simplified expression is only correct within a factor of three 
or four i. Since the efficiency of SFP depends on the product of neutrino magnetic moment 



^Assuming that condition (^Tj) is satisfied. Otherwise, the description of the LOW case is similar to that 
of LMA, and the accuracy of the approximate expression is as good as it is in the LMA case. 
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and magnetic field strength, only this product and not /i and B± separately can be probed 
by studying the solar z/g flux $p^. Comparing eqs. (^3]), and ( PUD we flnd that the 
sensitivity of the z/g flux to the product fiB± is strongest in the case of the VO solution of 
the solar neutrino problem; it is weaker in the case of the LOW solution and weakest for 
the LMA solution. This, however, does not necessarily mean that the LMA solution has 
the lowest sensitivity to the neutrino magnetic moment: the UeR flux in that case depends 
on the magnetic fleld in the core of the sun which may well be much stronger than the fleld 
in the convective zone, relevant for the VO case. 

We shall now discuss the present experimental upper bounds on as well as the 
sensitivity of the future experiments, and their implications. Currently, the most stringent 
upper bounds on come from the LSD experiment, < (1.7 x 10^^) $8^ at 90% CL [^, 
and from the Super-Kamiokande experiment, < (1.2 — 1.6) x 10^^ $8^ at 90% CL 



The Super-Kamiokande bounds were presented for several energy bins with E > 8 MeV. 
Future experiments are expected to improve these bounds (or discover the flux of solar UeR)'- 
KamLAND will be able to put a limit of 10~^ $8^ at 95% CL on the solar z/gR flux after one 
year of operation |5T|, and Borexino should be able to reach a similar sensitivity after a few 
years of data taking [^. The current limit ^ 1.5% $8^ implies, in the case of the LMA 
solution, a bound 

/i B, 



< w 



(37) 



_10-i2/i5 10 kC 

where B±{ri) is the average solar magnetic fleld in the neutrino production region, r < O.IRq. 
An experiment with the tritium radioactive antineutrino source has been recently proposed 
with the goal of putting an upper limit of 3 x lO^^^yU^ on the neutrino magnetic moment 
fi or measuring it if it is above this value ^3|]; if /i ~ 3 x lO'^^fi^ is found, the bound ( p7| ) 
would imply B±{ri) < 3 x lO''' G. Note that this limit is more stringent than the astrophysical 
one obtained from the requirement that the pressure of the solar magnetic fleld should not 
exceed the matter pressure {B± < 10^) ^4 . 



Conversely, if a rehable quantitative model of the solar magnetic fleld is developed, 
eq. (^) will limit the neutrino magnetic moment. For Bj_{ri) close to the above-mentioned 
astrophysical bound, the limit on /x would be /x < 10~^^ fiB, which is more than an order of 
magnitude more stringent than the expected limit from the planned laboratory experiment 



|33|. Unfortunately, no compelling model of the solar magnetic fleld exists at present. 

Similar considerations apply to the LOW and VO cases. For LOW, assuming that 
condition (pT]) is satisfled, the current limits on the solar z/g/j flux lead to 



Bn 



10-i2/i5 10 kC 



< 10 



parameterize the derivative of the solar magnetic fleld at 

-12, 



(3^ 



Rq, see 



where k and Bq 

eq. (|2|). For /i ~ 3 X 10"^^^, eq. limits this derivative to be \B'^{r)\R^ < 2.2 x 10^ 
kG/-R0. Conversely, if the actual value of |-Bj_(t')|_R0 is close to this value, the limit on /i 
from (p8D would be competitive with the expected upper bound from the planned laboratory 
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experiment [Q. If condition (|31|) is not satisfied, the preceding discussion of the LMA case 
apphes to the LOW welL 

In the VO case, the current upper bounds on imply 



/i B 
10^12^5 10 kG 



<5, (39) 



where B is the average magnetic field in the interval O.817i?0 < r < Rq defined in ([351) . 
For /X = 3 X 10~^2 this gives B < 17 kG. If, alternatively, some model considerations 
establish that the average field B is, for example, 100 kG, eq. (^) would lead to the limit 
< 5 X 10~^^. 

With the expected upper bound on the flux of solar z/gR from KamLAND, all the limits 
that we discussed above (eqs. (|37D - (0)) will be strengthened by about a factor of four. 

In the VO case, the flux of solar UeRS will have seasonal variations, similar to those of 
the z/eL flux. In the LOW and VO cases, for which the production is mainly driven by 
the magnetic field in the convective zone, the solar z/gR flux may also vary with time due to 
the 11-year variations of this magnetic field. For all the solutions, the solar Uf.R flux should, 
of course, also exhibit ~ 7% variations due to the variations of the distance between the 
sun and the earth. Low statistics may, however, make these variations difficult to detect. 
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Figure 2: Probabilities P{veL i^eR) corresponding to the UeL ^t^R amplitudes of eqs. (|l^) 
(solid curve) and ( p8| ) (dashed curve) for the LOW solution. Magnetic field decreases linearly from 
Bo = 10^ G at r = O.OSi?© to zero at r = Rq, fi = IQ-^Vb- 
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